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The  structure  and  propagation  of  oceanic  Yanai  waves  (also  known  as  mixed 
Rossby-gravity  waves)  are  investigated  by  analysis  of  satellite-derived  sea-surface- 
height  (SSH)  data.  Significant  spectral  peaks  along  the  dispersion  curves  of  first  and 
second  baroclinic  mode  Yanai  waves  are  identified  by  the  wavenumber -frequency 
spectral  analysis  of  SSH  for  the  period  10-17  days.  The  spatial  structure  and  prop- 
agation  of  these  modes  are  described  by  an  analysis  of  SSH  time  series  filtered  in 
the  frequency-wavenumber  domain  that  includes  strong  signals  along  the  disper¬ 
sion  curves  but  excludes  those  of  17-day  tropical  instability  waves.  The  difference  in 
meridional  structure  between  the  first  and  second  baroclinic  modes,  which  is  con¬ 
sistent  with  theory,  is  evident.  The  slow  eastward  propagation  of  first  and  second 
baroclinic  mode  Yanai  wave  packets,  which  is  consistent  with  the  group  velocity 
derived  from  theory,  is  also  isolated  from  the  SSH  data.  Year-to-year  variability 
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1.  Introduction 

Yanai  waves  (often  referred  to  as  mixed  Rossby-gravity 
waves)  play  an  important  role  in  driving  a  variety  of 
phenomena  in  the  tropical  ocean  and  atmosphere  (e.g.  Yanai 
and  Hayashi,  1969;  Zangvil,  1975;  Liebmann  and  Hendon, 
1990;  Takayabu  and  Nitta,  1993;  Lyman  etaL>  2007).  The 
structure  and  dispersion  relationships  of  Yanai  waves  were 
derived  as  one  of  the  eigen -modes  of  the  linear  shallow- 
water  equation  on  the  equatorial  ^-plane  in  the  1960s 
(Matsuno,  1966).  The  wave  exhibits  an  eastward  group 
velocity  at  all  wavenumbers.  It  behaves  like  a  westward 
propagating  Rossby  wave  in  the  large  negative  wavenumber, 
whereas  it  has  a  structure  similar  to  an  inertia  gravity  wave 
in  the  large  positive  wavenumber.  Hence,  for  the  negative 
wavenumber,  the  wave  crest  and  trough  propagate  westward 

(£>  Published  in  201 1  by  John  Wiley  &  Sons  Ltd. 


whereas  the  wave  energy  (packet)  propagates  eastward  with 
the  speed  of  group  velocity.  Observational  evidence  of  Yanai 
waves  was  presented  by  the  analysis  of  winds  in  the  lower 
stratosphere  almost  at  the  same  time  as  its  theoretical 
derivation  (Yanai  and  Maruyama,  1966).  Since  then,  many 
studies  discussed  its  structure,  propagation,  and  generation 
and  forcing  mechanisms  for  both  the  atmosphere  (e.g.  Itoh 
and  Ghil,  1988;  Hendon  and  Liebmann,  1991;  Takayabu, 
1994)  and  ocean  (e.g.  Wunsch  and  Gill,  1976;  Moore  and 
Philander,  1977;  Weisberg  et  al>  1979;  Zhu  etal ,  1998). 

Using  outgoing  long- wave  radiation  (OLR)  data,  eastward 
group-propagation  of  the  atmospheric  Yanai  waves  coupled 
to  convection  was  identified  by  the  combination  of 
wavenumber- frequency  spectral  analysis  and  regression 
analysis  (Wheeler  and  Kiladis,  1999;  Wheeler  2000). 
The  period  of  typical  atmospheric  Yanai  waves  is  2-6  days, 
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and  the  group  velocity  is  about  5ms*"1,  which  is  shown  to 
he  well  resolved  by  the  twice-daily  OLR  data  (Wheeler  and 
Kiladis,  1999;  Wheeler  etaiy  2000).  Although  the  period 
of  oceanic  Yanai  wave  is  longer  (about  10-20  days),  the 
observational  evidence  of  its  eastward  group-velocity  has  not 
been  reported  because  until  recently  the  spatial  and  temporal 
resolutions  of  satellite  and  in  situ  data  were  insufficient 
to  explicitly  describe  the  propagation  of  the  wave  packet. 
This  situation  has  begun  to  change  because  of  the  multiple 
satellite  launches  in  recent  years,  and  the  most  recent  satellite 
altimetry  products  are  shown  to  he  able  to  well  resolve  the 
signal  of  oceanic  equatorially  trapped  waves  including  Yanai 
waves  and  to  describe  their  observed  dispersion  relations 
(Wakata,  2007;  Farrar,  2008,  2011;  Shinoda  etaly  2009; 
Shinoda,  2010). 

Using  the  sea-surface  height  (SSH)  data  from  recent 
satellite  altimeter  measurements,  Shinoda  (2010)  identified 
prominent  spectral  peaks  along  the  Yanai  wave  dispersion 
curves  for  the  first  and  second  baroclinic  modes  by 
wavenumber -frequency  spectral  analysis.  The  spatial 
structure  and  the  phase  relationship  of  17-day  tropical 
instability  waves  (TIWs;  Lyman  et  aly  2007)  are  described  by 
the  space  and  time  filtered  SSH  and  sea-surface  temperature 
(SST),  showing  that  the  structure  of  the  17-day  TIWs  is 
similar  to  that  of  Yanai  waves  derived  from  theory.  The 
meridional  structure  of  SSH  fluctuation  associated  with 
17-day  TIWs  is  shown  to  be  anti -symmetric  with  regard 
to  the  Equator.  The  similar  anti-symmetric  structure  is 
also  found  in  the  SST  variation  associated  with  17-day 
TIWs.  However,  his  analysis  did  not  clearly  demonstrate 
the  difference  in  structure  between  the  first  and  second 
baroclinic  modes  because  it  focused  on  showing  the  structure 
of  TIWs  for  the  period  around  17  days,  where  the  spectral 
peaks  along  the  two  dispersion  curves  are  merged  in  the 
frequency-wavenumher  space. 

In  this  study,  we  investigate  the  spectral  characteristics, 
structure  and  propagation  of  Yanai  waves  based  on  the 
combination  of  wavenumber -frequency  spectral  analysis 
and  cross-correlation  (regression)  analysis  using  the 
most  recent  version  of  SSH  product  from  satellite 
altimeter  measurements.  The  SSH  data  are  filtered  in 
the  frequency- wavenumber  domain  where  first  and 
second  mode  signals  are  clearly  distinguished.  Unlike  the 
atmospheric  Yanai  wave,  the  signal  of  the  oceanic  Yanai 
wave  can  be  clearly  separated  into  its  distinct  vertical  modes, 
and  thus  we  are  ahle  to  descrihe  the  spatial  distribution  and 
structure  of  first  and  second  baroclinic  mode  Yanai  waves 
separately.  In  addition,  the  slow  eastward  propagation  of 
first  and  second  mode  Yanai  wave  packets  are  isolated, 
and  the  propagation  speeds  of  these  modes  are  shown 
to  be  consistent  with  the  group  velocities  obtained  from 
theoretical  dispersion  curves. 

2.  Wavenumber -frequency  spectral  analysis 

We  use  daily  multisatellite  analysis  of  SSH  anomalies  on  a 
1  x  1°  grid  over  the  period  2003-2009.  These  analyses  are 
obtained  from  Archiving,  Validation  and  Interpretation  of 
Satellite  Oceanographic  (AVISO)  data.  It  should  be  noted 
that  the  1°  x  1°  gridded  daily  SSH  products  are  a  smoothed 
version  of  the  measurements,  which  are  derived  by  objective 
analysis,  and  that  the  actual  resolution  of  the  data  is  still 
uncertain  (e.g.  Chelton  and  Schlax,  2003).  Although  the 
daily  1°  x  1°  SSH  products  are  currently  available  from 


October  1992,  we  use  only  the  recent  data  from  2003  when 
measurements  from  multiple  satellites  were  included  in  the 
ohjective  analysis. 

Previous  studies  have  demonstrated  that  the  space  and 
time  structure  of  equatorially  trapped  waves,  including  Yanai 
waves,  are  effectively  isol  ated  by  the  wavenu  mber-  frequency 
spectral  analysis  for  both  the  atmosphere  and  ocean.  Here, 
we  first  repeat  the  calculation  of  wavenu  mber- frequency 
spectra  performed  by  Shinoda  (2010)  using  the  most  recent 
version  of  AVISO  SSH  product.  Since  the  structure  of 
the  Yanai  wave  is  anti-symmetric  about  the  Equator,  the 
anomalies  of  SSH  are  first  decomposed  to  anti-symmetric 
and  symmetric  components  (Wheeler  and  Kiladis,  1999). 
Then  a  complex  fast  Fourier  transform  (FFT)  of  SSH 
anomalies  is  performed  over  longitude  to  decompose  SSH 
into  zonal  wavenumbers,  and  an  additional  FFT  is  used  to 
determine  the  spectral  power  of  each  resolved  wavenumber 
by  frequency.  Further  details  of  the  analysis  are  found  in 
Wheeler  and  Kiladis  (1999)  and  Shinoda  et  at.  (2009). 

Figure  1  shows  contours  of  the  base- 1 0  logarithm  of  power 
in  wavenu  mber- frequency  space  calculated  from  the  anti¬ 
symmetric  component  of  SSH  anomalies  in  the  equatorial 
area  (5°N-5°S)  in  the  Pacific  Ocean  (east  of  1503E).  Solid 
lines  in  Figure  1  are  the  dispersion  relationship  of  Yanai 
waves  for  the  first  and  second  baroclinic  modes: 

k  =  co/c  —  ft/co  (1) 

where  co  is  the  angular  frequency,  k  is  the  wavenumher,  c 
is  the  internal  gravity  wave  speed,  and  ft  is  the  variation 
of  the  Coriolis  parameter  with  latitude.  Following  Wheeler 
and  Kiladis  ( 1999)  and  Shinoda  et  a/. (2009),  the  background 
spectrum  is  removed  from  the  original  spectrum.  First,  the 
raw  spectra  of  positive  and  negative  wavenumber  domains 
are  averaged  to  make  the  spectrum  symmetric  with  regard  to 
wavenumber  0.  Then  the  background  spectrum  is  calculated 
by  applying  a  1-2-1  filter  multiple  times  in  frequency  and 
wavenumber.  Here,  the  number  of  passes  of  the  12-1  filter 
is  seven  in  frequency  and  four  in  wavenumber. 

Prominent  spectral  signals  of  Yanai  waves  are  evident, 
with  their  peak  corresponding  to  the  dispersion  curve  of  an 
equivalent  depth  0.8  m.  This  equivalent  depth  corresponds 
to  equatorial  Kelvin  waves  with  a  phase  speed  of  ~2.8  ms-1, 
which  is  consistent  with  the  first  baroclinic  mode  Kelvin 
waves  in  the  equatorial  Pacific  Ocean  (e.g.  Johnson  and 
McPhaden,  1 993;  Picau  t  and  Sombardier,  1 993;  Kessler  et  ai , 
1995;  Cravatte  eta/.,  2003;  Shinoda  et  a/.,  2008).  Signals  of 
the  second  baroclinic  mode  Yanai  wave  are  also  found  for 
the  equivalent  depth  of  0.26  m,  which  is  consistent  with 
the  second  mode  Kelvin  wave  phase  speed  of  M.6  m  s“l 
(e.g.  Cravatte  et  a/.,  2003).  The  spectral  signals  of  the  second 
baroclinic  mode  are  comparahle  to  those  for  the  first  mode. 

In  order  to  describe  the  spatial  variation  of  Yanai  wave 
activity,  space  and  time  filtering  (via  spectral  transform) 
centred  around  the  spectral  peak  of  each  mode  (the  area 
of  the  two  boxes  in  Figure  1)  is  applied  to  the  SSH  time 
series.  These  areas  for  the  filter  are  chosen  ro  clearly  separate 
the  first  and  second  baroclinic  modes.  Although  significant 
peaks  are  found  around  the  period  of  17  days  because 
of  tropical  instability  waves  (Lyman  et  a/.,  2007;  Shinoda, 
2010),  first  and  second  mode  dispersion  curves  are  merged 
around  this  period.  Hence  most  areas  around  this  period  are 
excluded  by  the  filtering.  Note  that  the  time  series  are  not 
decomposed  to  symmetric  and  anti-symmetric  components 
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Figure  I.  Zonal  wavenumber -frequency  power  spectrum  of  sea-surface- 
height  anomalies.  The  colour-bar  number  indicates  the  base- 10  logarithm. 
The  background  power  is  removed.  The  curves  indicate  Yanai  wave 
dispersion  relations  for  equivalent  depths  of  0.8  m  and  0.26  m. 


for  the  filtering.  Figure  2(a)  shows  the  standard  deviation 
of  SSH  time  series  at  each  location,  which  are  filtered 
in  space  and  time  for  the  area  surrounded  by  the  upper 
box  in  Figure  1.  This  area  covers  the  significant  spectral 
peak  associated  with  first  baroclinic  mode  Yanai  waves. 
The  maximum  of  SSH  fluctuation  is  found  around  3'S, 
between  the  dateline  and  140UW,  and  there  is  another  weak 
local  maximum  in  the  Northern  Hemisphere  around  2°N, 
between  1 50°  W  - 1 20°  W. 

Figure  2(b)  shows  the  standard  deviation  of  SSH  time 
series  for  the  second  baroclinic  mode,  which  are  filtered  in 
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Figure  2.  (a)  Standard  deviation  of  sea-surface-height  anomalies  (cm) 
filtered  at  periods  and  wavelengths  in  the  area  surrounded  by  the  upper 
box  shown  in  Figure  1 .  (b)  Same  as  (a)  except  for  the  lower  box  in  Figure  1 . 


space  and  time  for  the  area  indicated  in  the  lower  box  in 
Figure  1.  The  spatial  pattern  of  SSH  fluctuation  is  similar  but 
notably  different  from  that  for  the  first  baroclinic  mode.  The 


SSH  Lag  0  Mode  1 


180  175W  I70W  165W  I60W  I55W  150W  H5W  HOW  135W  130W  125W  120W' 

SSH  Lag  0  Mode  2 


0  16 
|  0.12 
o.oe 

0.04 

0 

-0.04 

I  -0.08 


k  0.2 
0.16 
0. 1 2 
0.06 
0.04 
0 

“0  04 
-0  08 
-0  12 


Figure  3.  Upper  panel:  sea-surface-height  (SSH)  anomalies  regressed  onlo space-  and  lime-HIlered  SSH  anomalies  al  2.5  S,  15Q°W  for  the  first  baroclinic 
mode  Yanai  wave.  Lower  panel:  same  as  the  upper  panel  except  for  the  second  baroclinic  mode. 
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maximum  amplitude  is  found  around  2°S,  180  -120°W. 
The  weaker  local  maximum  in  the  Northern  Hemisphere  is 
also  found  around  1CN-2°N,  between  150°W— 120°W. 

These  spatial  distributions  of  Yanai  wave  activity  are 
quite  different  from  those  associated  with  17 -day  tropical 
instability  waves  (TlWs),  in  which  the  location  of  the 
maximum  activity  is  centred  around  120°W  in  the  eastern 
Pacific  (Shinoda,  2010).  Because  the  minimum  of  mean  SST 
is  found  on  the  Equator  in  the  eastern  Pacific  and  thus  the 
sign  of  meridional  SST  gradient  is  opposite  in  the  Northern 
and  Southern  Hemispheres,  the  meridional  component  of 
velocity  fluctuation  associated  with  Yanai  waves  generates 
the  SST  variation  that  is  anti-symmetric  with  regard  to  the 
Equator.  Therefore,  significant  spectral  signals  of  SST  are 
evident  around  the  17-day  period  along  the  Yanai  wave 
dispersion  curves.  However,  the  meridional  SST  gradient  is 
much  weaker  in  the  central  Pacific  where  strong  activity  of 
Yanai  waves  with  a  10- 16-day  period  is  observed,  and  thus 
significant  spectral  peaks  along  the  Yanai  wave  dispersion 
curves  for  these  periods  are  not  found  in  the  SST  spectrum 
(Shinoda,  2010). 

The  result  indicates  that  the  meridional  length  scale 
for  the  first  baroclinic  mode  is  larger  than  that  of  the 
second  mode,  which  is  at  least  qualitatively  consistent  with 
theory.  This  difference  in  meridional  structure  was  not 
clearly  demonstrated  for  Yanai  waves  with  a  17-day  period 
(Shinoda,  2010).  As  the  two  dispersion  curves  around  this 
period  (17  days)  are  very  close  in  frequency-wavenumber 
space  (see  the  two  boxes  of  figure  1(b)  in  Shinoda,  2010), 
signals  of  two  modes  may  not  be  well  isolated  by  the  analysis 
of  Shinoda  (2010). 

The  observed  meridional  structure  shows  a  larger 
amplitude  in  the  Southern  Hemisphere,  which  is  not 
explained  by  a  linear  wave  solutions  derived  from  shallow- 
water  equations.  The  existence  of  strong  mean  currents  in 
this  region  may  influence  the  structure  (Lyman  et  al.y  2007). 
Also,  the  spatial  resolution  of  the  SSH  data  is  not  sufficient 
for  further  quantitative  comparison  with  theory. 

3.  Lagged  regression  analysis 

In  order  to  further  examine  the  property  of  Yanai  waves, 
such  as  the  structure  and  propagation  of  individual  waves, 
time  series  of  unfiltered  SSH  anomalies  are  regressed  against 
space-  and  time-filtered  SSH  anomalies  at  2.5~"S,  150°W 
where  the  large  SSH  variance  is  found  in  both  first  and 
second  baroclinic  modes  (Figure  2(a)  and  (b)).  Note  that  the 
data  are  not  decomposed  to  symmetric  and  anti-symmetric 
components  for  the  regression  analysis.  The  spatial  pattern 
of  SSH  fluctuation  (Figure  3(a)  and  (b))  for  both  first 
and  second  baroclinic  modes  is  nearly  anti-symmetric 
with  regard  to  the  Equator,  which  is  consistent  with  the 
meridional  structure  derived  from  theory. 

Figure  4  shows  the  time-longitude  (Hovmoller)  diagram 
of  the  regressed  SSH  along  3°S  for  the  first  and  second 
baroclinic  modes.  While  the  westward  propagation  of 
individual  Yanai  waves  is  evident  for  both  modes,  the  packet 
of  these  waves  slowly  propagates  eastward.  The  speed  of  the 
eastward  propagation  for  the  first  baroclinic  mode  is  faster 
than  the  second  baroclinic  mode.  These  signals  of  slow 
eastward  propagation  are  consistent  with  the  group  velocity 
of  Yanai  waves  for  the  first  baroclinic  mode  (1.03  m  s-*1) 
and  second  baroclinic  mode  (0.64  ms”1)  derived  from  the 
dispersion  curves  in  Figure  1.  Since  unfiltered  SSH  fields  are 


Figure  4.  Upper  panel:  longitude-lime  diagram  of  regressed  sea  surface- 
height  anomalies  for  the  first  baroclinic  mode.  The  solid  line  indicates  the 
group  velocity  (1.03  m  s"1)  for  the  first  baroclinic  mode  Yanai  wave  at  ihc 
wavenumber  4007  km  (ihe  global  wave  number  —  1 0  for  the  horizontal  axis 
in  Figure  I)  derived  from  the  dispersion  curve  in  Figure  1.  lower  panel: 
same  as  the  upper  panel  except  for  the  second  baroclinic  mode  of  which 
lhe  group  velocity  is  0.64  m  s”1 . 


regressed  against  the  Yanai  wave  filtered  time  series  at  the 
base  point,  only  a  real  associated  signal  that  is  coherent  with 
the  base  point  time  series  is  extracted.  Hence,  the  analysis 
demonstrates  that  signals  of  slow  eastward  propagation  of 
the  Yanai  wave  packet  can  be  detected  from  the  satellite- 
derived  SSH  data. 

Longitude-time  diagrams  of  space- time  filtered  SSH 
are  shown  in  Figure  5  in  order  to  examine  the  temporal 
variation  of  the  Yanai  wave  amplitude  for  each  vertical  mode. 
While  the  amplitudes  of  first  and  second  baroclinic  modes 
are  comparable  during  this  period  (mid-May  2004  to  July 
2005),  the  relative  amplitude  of  each  mode  depends  on  the 
period.  The  first  mode  Yanai  waves  have  larger  amplitude 
during  June-September  2004,  whereas  the  amplitude  of 
the  second  mode  is  larger  during  January- June  2005.  It 
should  he  noted  that  the  filtered  time  series  for  the  same 
frequency-wavenumber  domain  display  similar  eastward 
propagation  of  the  wave  envelope  for  any  field  (e.g.  red- 
noise  fields)  because  of  the  design  of  the  filter,  but  the 
temporal  variation  of  wave  amplitude  for  each  baroclinic 
mode  can  be  identified  by  these  filtered  time  series. 

4.  Year-to-year  variation 

Shinoda  etui  (2009)  indicated  that  the  interannual  variation 
of  Rossby  wave  and  tropical  instability  wave  activity  in 
the  eastern  equatorial  Pacific  is  moderately  correlated  with 
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Figure  5.  Left  panel:  longitude-time  diagram  of  sea-surfacc-height  anomalies  filtered  in  the  frequency-wavenumber  domain  indicated  in  the  upper  box 
in  Figure  I  for  the  period  14  May  2004  to  7  August  2005.  The  solid  straight  line  indicates  the  group  velocity  (1.03  m  )  at  the  wavenumber  4007  km. 
Right  panel:  same  as  the  left  panel  except  that  the  domain  of  the  filter  is  indicated  in  the  lower  box  in  Figure  1. 


an  El  Nino-Southern  Oscillation  (ENSO)  index  based  on 
the  analysis  of  SSH  and  SST  time  series  filtered  in  the 
frequency- wavenumber  domain.  Here,  the  same  method 
as  used  in  Shinoda  et  al.  (2009)  is  applied  to  investigate  the 
year-to-year  variation  of  Yanai  wave  activity.  The  intensity 
of  Yanai  wave  activity  is  measured  as  the  standard  deviation 
of  the  space-  and  time- filtered  SSH  time  series  displayed  in 
Figure  5.  Figure  6  shows  the  monthly  time  series  of  first 
and  second  baroclinic  mode  Yanai  wave  activity  for  the 
area  1 80J- 1 10C  W,  4°S-0°,  where  a  significant  amplitude  of 
Yanai  wave  is  observed.  Time  series  of  first  and  second  modes 
have  some  similarities  and  also  notable  differences,  with  a 
moderate  correlation  between  the  two  modes  (correlation 
coefficient:  i?  =  0.54).  For  example,  activity  for  both  first 
and  second  modes  is  strong  during  late  2004  and  relatively 
weak  during  early  2009.  Also,  a  rapid  increase  of  the  activity 
is  evident  for  both  first  and  second  modes  during  2009. 
However,  the  intensity  of  the  second  mode  is  about  twice  as 
large  as  the  first  mode  during  most  of  the  periods  in  2008. 


Further  quantitative  comparisons  between  first  and  second 
modes  may  not  be  appropriate  since  the  altimeter  product 
might  better  resolve  the  signal  of  a  second  mode  with  a  wave 
period  that  is  longer  than  the  first  mode. 

The  time  series  of  Yanai  wave  activity  are  then  compared 
with  ENSO  indices.  Since  it  has  become  widely  accepted 
that  there  are  two  different  types  of  El  Nifto  in  recent 
years  (Larkin  and  Harrison,  2005;  Ashok  etal ,  2007;  Kao 
and  Yu,  2009;  Kug  etal,  2009),  two  indices  are  used  for  the 
comparison.  These  are  the  Nino3  ( 1 50^  W-90°  W,  5°N-  5US ) 
SST  anomaly  and  the  ENSO  Modoki  Index  (EMI;  Ashok 
et  al. ,  2007)  defined  as 

EMI  =  [SST]a  -  0.5  x  [SST]b  -  0.5  x  [SST]C  (2) 

where  [SST] A,  [SST]B  and  [SST]C  axe  the  average  SST 
anomalies  over  regions  A  (165°E-140°W,  10°S-10°N), 
B  (110CW-70'W,  15JS-5CN)  and  C  (125°E- 145°E, 
10°S-20’N)  respectively.  Monthly  SST  data  from  the 
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Figure  6.  (a)  The  monthly  mean  intensity  of  sea  surface- height  (SSH) 
variability  associated  with  first  (solid  line)  and  second  (dashed  line) 
baroclinic  mode  Yanai  wave  for  the  area  180  -110’ W,  4°S-0°.  The 
intensity  is  measured  by  the  standard  deviation  of  space-  and  lime-filtered 
SSH  time  series,  (b)  Summation  of  the  two  curves  in  (a)  (solid  line), 
ENSO  Modoki  Index  (EMI;  dashed  line:  see  text  for  the  detail),  and  Nino  3 
sea-surface  temperature  anomaly  (doited  line).  Time  series  are  normalized 
by  their  standard  deviation  and  the  means  are  subtracted.  This  figure  is 
available  in  colour  online  at  wilcyonlinelibrary.com/journal/qj 


blended  product  of  satellite/m  sifu  observations  (Reynolds 
a  al.y  2002)  are  used  to  calculate  EMI  and  Nino3  SST 
anomalies.  The  new  type  of  El  Nino  is  often  referred  to 
as  El  Nifio  Modoki  or  Central  Pacific  El  Nino  where  the 
maximum  SST  anomaly  is  located  in  the  central  Pacific,  and 
is  Hanked  by  anomalous  cold  SSTs  at  the  western  and  eastern 
ends  of  the  equatorial  Pacific. 

Figure  6(b)  compares  the  Yanai  wave  activity  with  the 
Nino3  SST  anomaly  and  EMI.  While  the  time  series  of 
Yanai  wave  activity  is  not  similar  to  the  index  of  the 
traditional  El  Nino  ( R  —  0.25),  it  is  significantly  correlated 
with  EMI  ( R  =0.65).  The  EMI  is  better  correlated  with 
the  first  baroclinic  mode  Yanai  wave  activity  ( R  =  0.68) 
than  the  second  mode  ( R  =  0.40).  Whereas  the  maximum 
anomalies  of  SST  are  located  in  the  eastern  Pacific  during 
the  traditional  El  Nino,  major  SST  anomalies  during  El  Nino 
Modoki  events  are  located  in  the  central  Pacific  where  Yanai 
wave  activity  is  relatively  strong.  Hence,  large-scale  wind 
anomalies  associated  with  anomalous  SSTs  in  the  central 
Pacific  may  generate  anomalous  ocean  currents  different 
from  those  during  the  traditional  El  Nino. 

A  recent  study  by  Shinoda  el  al.  (2011)  demonstrates  that 
anomalous  zonal  currents  near  the  Equator  show  conver¬ 
gence  (divergence)  with  eastward  (westward)  currents  west 
of  the  dateline  and  westward  (eastward)  currents  east  of  it 


during  the  2004  (2008)  El  Nino  (La  Nifia)  Modoki  event. 
Thus  the  westward  South  Equatorial  Current  (SEC)  east 
of  the  dateline  is  enhanced  (reduced)  during  2004  (2008) 
El  Nino  (La  Nina)  Modoki.  The  stronger  SEC  (and  thus 
stronger  horizontal  shear)  during  2004  El  Nino  Modoki 
might  have  an  impact  upon  the  subseasonal  variance  of 
meridional  velocity  (and  thus  generation  of  Yanai  waves) 
through  enhanced  instability  of  equatorial  currents.  Also, 
subseasonal  (5-20  day  period)  wind  activity  is  enhanced 
in  the  central  Pacific  east  of  the  dateline  during  the  2004 
El  Nifio  Modoki  event  in  association  with  the  warm  SST 
anomaly  (not  shown).  Subseasonal  winds  could  affect  Yanai 
wave  activity  through  directly  generating  the  meridional 
component  of  surface  currents.  However,  the  time  series 
discussed  here  are  not  long  enough  to  investigate  the  inter¬ 
annual  variation  associated  with  ENSO  in  greater  detail. 
Further  observational  and  modelling  studies  are  necessary 
to  fully  understand  the  physical  processes  responsible  for 
the  significant  correlation  between  Yanai  wave  activity  and 
the  ENSO  Modoki. 

5.  Discussion  and  conclusions 

Unlike  the  atmospheric  equatorial  waves,  signals  of  oceanic 
waves  are  clearly  separated  into  distinct  vertical  modes. 
Therefore  the  observed  structure  and  group  velocity  for 
first  and  second  baroclinic  mode  oceanic  Yanai  waves  can 
be  compared  with  theory  separately,  and  the  differences 
between  the  two  modes  can  be  described  explicitly.  This  is 
not  possible  for  the  atmospheric  Yanai  wave  as  the  distinct 
equivalent  depths  of  different  vertical  modes  cannot  be 
clearly  identified  because  of  the  absence  of  a  solid  boundary 
at  the  top  of  the  troposphere  and  the  interaction  between  the 
waves  and  moist  convection.  Accordingly,  the  broad  cover¬ 
age  of  recent  satellite  SSH  measurements  provides  excellent 
data  sets  to  investigate  the  characteristics  of  observed 
equatorially  trapped  waves  such  as  the  phase  and  group 
velocities  of  Yanai  waves  and  their  comparison  with  theory. 

Because  the  time  and  space  resolutions  of  recent  satellite 
altimetry  products  are  substantially  improved,  it  is  now 
feasible  to  describe  space  and  time  structure  of  oceanic  Yanai 
waves,  including  the  identification  of  the  spatial  distribution 
of  wave  activity  and  the  depiction  of  group  velocity.  Our 
statistical  analyses  clearly  demonstrate  that  signals  of  Yanai 
waves  can  now  be  detected  in  daily  SSH  products  derived 
from  the  recent  satellite  measurements.  First  and  second 
baroclinic  mode  Yanai  waves  can  be  clearly  distinguished, 
and  their  spatial  structure  can  be  isolated.  Also,  the  group 
velocity  for  each  mode  can  be  explicitly  described  as  a 
slow  eastward  propagation  of  the  wave  packet.  Year-to-year 
variability  of  Yanai  wave  activity  is  compared  with  that  of 
the  ENSO  based  on  the  analysis  of  space  and  time- filtered 
SSH  time  series,  showing  a  moderate  correlation  with  the 
ENSO  Modoki.  These  statistical  representations  of  oceanic 
Yanai  waves  are  important  for  understanding  their  dynamics 
and  generation  mechanism,  and  also  useful  for  evaluating 
numerical  model  performances. 
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